DEMIMARTINGALES AND THE FUNCTIONAL HILL PROCESS 
FOR SMALL PARAMETERS 

ADJA MBARKA FALL, GANE SAMB LO, AND CHEIKHNA HAMALLAH NDIAYE 

Abstract. Association of random variables and Demimartingales are recent 
fields for handling asymptotic behaviors of sums of dependent random vari- 
ables. We apply their techniques to establish the asymptotic law of the demi- 
martingale 

fc-i fe-i fc-i 

^/(i)(cxp(- 7 J2 E h /h)-cxp(-yJ2 E h/h)), 
3=1 h=j+l h=j 

where E\, E2, ... arc indcpcndant standard exponential random variables, 7 > 
0, k is a positive integer and f(j) is an increasing function of the integer j > 1. 
We next apply the results to find the asymptotic behavior the functional Hill 
process for small parameters whithin the Extreme Value Theory (EVT) field. 
Such a result would have been very hard to find whithout demimartingales 
techniques. 



I. Introduction 

We are concerned in this paper with some demimartingale and its 
application in extreme value theory. Namely let n > 1, 1 < k(n) < n, 
7 > 0, let Ei tn , E 2t n, E n , n be a sequence of independent standard 
exponential random variables defined on the same probability space 
(fi, A, P), let f(j) be a nondecreasing function of j e M\{0}, and finally 
let 

fc(n)-l 

(1.1) W k{n) M) = E /0')(exp(F i+1>B ) - exp(F j , n )), 

j'=i 

where 

fe-i 

Fj, n = -771%/^ ^fc,™ = 0. 
h=j 

In this note, we intend to use the relatively recent theory of associated 
random variables and demimartingales to determine the limiting law 
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W k (n),n(f) when 

fc-i 

(ki) ^(/(j) - f(j - i))r 1/2 < oo, 

and next to apply the results to solve an open problem in statistics of 
the extreme value Theory. We then may begin to introduce to associ- 
ated random variables which goes back to Lehmann (1966) [11] in the 
bivariate case. 

The concept of associated variable generalizes that of independence 
and seems to model a great variety of stochastic models. We point 
out now that we will equally say that a sequence of random variables 
is associated or that the elements of the sequence are associated. As 
immediate examples, Gaussian random vectors with nonnegatively cor- 
related components(see Pitt [12] ). Also the order statistics associated 
with a finte set of independent are associated . As well a homogenuous 
Markov chain is also associated (Daley [2]). Such a property also arises 
in Physics, and is quoted under the name of FKG property (Fortuin, 
Kastelyn et Ginibre (1971) [7]), in percolation theory and even in Fi- 
nance (see Pan Jiazhu [16J). The final definition is given by Esary, 
Proschan et Walkup (1967) [6] in 

Definition 1. A finite sequence ofrv's (Xi, ...,X n ) are associated when 
for couple of real and coordinate-wise nondecreasing functions h et g 
defined on M. n , we have 

(1.2) Cov(h(Xi, X n ), g(X 1 ,...,X n ))>0 

An infinite sequence of rv's are associated whenever all its finite sub- 
sequences are associated. 

We have a few number of interesting properties to be found in Rao 
([13]) : (PI) A sequence of independent rv's is associated. (P2) Par- 
tial sums of associated rv's are associated. (P3) Nondecreasing func- 
tions of associated variables are associated. (P4) Let the sequence 
Zi, Z 2 , Z n be associated and let (aj)i<j< n be positive numbers and 
(bi)i<i<n rea l numbers. Then the rv's ai{Zi — bi) are associated. 

Demimartingales are set from associated centered variables exactly 
as martingales are derived from partial sums of centered independent 
random variables. We have 
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Definition 2. A sequence of rv's {S n ,n > 1} in P) is a 

demimartingale when for any j > 1,/or any coordinatewise nondecreas- 
ing function g defined on IR d , we have 

(1.3) E((S j+1 -Sj) g(S 1 ,...,S j )) >0, j > 1. 

Two particular cases should be pointed out. First any martingale is a 
demimartingale. Secondly, partial sums S n = Xi + ...+X n of associated 
and centred randoms variables Xi,X 2 ,... form a demimartingale for, 



in this case, (1.3) becomes : 

E {(S j+1 - Sj) g(S 1 , Sj)} — E {X j+1 g^, Sj)} = Cov {X j+1 g(Sx, Sj)} , 

since EXj + i = 0. Since (xi, ...,Xj + i) i — )• a; J+ i et (xi, Xj+%) 1 — ^ 
g(xi, ...,Xj) are coordinate-wise nondecreasing functions and since the 
Xi,X 2) .. are associated, we get 

E {(S J+1 - 5,-) (/(Si, Sj)} = Cov {X j+1 g(S l , Sj)} > 0. 

We claim that this thoery is enough to handle the limiting law issue 
of Wk, n and its application in statistics of the extreme to solve open 
problems. We organize the remainder of the paper as follows. In Sub- 
section 2.1 of Section |2j we link Wk, n to demimartingales theory and 
provide our main result on Wk,n with the help of recent theorems on 
demimartingales in Subsection 22. In Section [3j we give some appli- 



cations on the functional Hill process and concluding remarks. The 
paper is finished by an appendix about moment computations. 

2. Main results and proofs 
2.1. How demimartingales may help? 

2.1.1. Preliminary lemma. 

Lemma 1. For each n > 1, the family {E(Wk t n) ~ Wk,m 1 < k < n} 
is a demimartingale. 

2.1.2. Proof. We have 

w Kn = f(k - 1) - (fl TO) - M - !)) ex P(-^ E E Pl h )) 

\ j=i h=j ' 

j=k-l 

= f(k - 1) - fU)Sik,n 

with s* An = eM-iEl^E^/h), f(j) = /O0-/(i-i)J(o) = o. 
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We have by (14.1. 11) 



fc-i 



*h = E ( S lk,n) = exp(- ^log(l + 1 /h)) 

h=j 



and 
(2.1) 

and next 
(2.2) 



E{W, 



k 



j=k-l 

f(k - 1) - 7UKk,n 
3=1 



w, 



k.n 



E(W t 



j i ) 



W, 



k,n 



J E /(;)(** -4) 

3=1 

which is a sum of centred associated rv's. By the properties (P1)-(P6) 
above, the f(j)(S* — s*) are associated since f(j) > for all j > 1. 
Then their partial sums form a demimartingale. 

2.2. The main result. 

Theorem 1. Let f(j) be an increasing function of the integer j > 1 
such that 

oo 

(2.3) »./)./ ' 2 - 'x. 

3=1 

and let 1 < k < n. Then Wk t n(f) converges p.s to a finite random 
variable W(f) such that E \W(f)\ < oo. Further if f(j) = j T , < r < 
1/2, then 

j=k-l 
3=1 

also converges to a random variable W{r) such that E \ W{r)\ < oo. 

2.3. Proof of Theorem [Tj We begin to consider an infinite sequence 
of independent and standard exponential random variables E\,E2,... 
defined on the same probability space. Under the hypotheses of the 
theorem, we have 

fc-i fe-i fc-i 

W:(f) = J2fti)(eM-"f E E h /h)-eM-lY. E h/h)),k>l, 

j=l h=j+l h=j 

with the convention Ylh^k-^h/h — 0. We have for any any function /, 

1 < k < n, 

W* k {f) = d WUf), 
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where =d stands for the equality in distribution. It is clear that if we 
prove that W k (f) converges a.s. to a finite random variable W(f) as 
k —¥ oo, we will able to conclude that W kjn (f) converges in distribution 
to W(f) as k = k(n) — » +00. Let us now prove that W k (f) converges 
a.s. to a finite random variable W(f) as k — > 00 by applying The- 
orm 2.4.2 of Rao (|13j) for demimartingales since W k (f) also forms a 
demimartingale. Is suffices to check that 

(2.4) lim sup E \W k \ < 00 

Denote S£ fc = exp(- 7 £*l} E h /h), 1 < j < fc - 1 and S* kjk = 1. By 
Cauchy-Schwarz's inequality and next by Minkoswki's one, we have 

fc-i 

E \W* k \ < {EWf) 1/2 = \\W* k \\ 2 < \\fUm k - s* fe )|| 2 

J=l 

fc-1 

<E^')(V r ar(^ fe )) 1 /2. 



Now by (4.2), fiNVariSU)) 1 ' 2 is dominated by Const. f(j)j~ 1/2 for 



large values of j. Thus Rao's condition (2.4)is valid and this gives the 
first part of the proof. Now for f(j) = j T , 

JU)r 1/2 ~ Const. r ((3/2)-) 

and j : ~(( 3 / 2 )~ T ) converges if and only if r < 1/2. This gives the 
second part of the Theorem [TJ 

2.4. Finite-dimensional limiting law. Consider the space J 7 of func- 



tions increasing functions f(j) of j > 1 satisfying (2.4) and consider 
{/1, f r } C F, r > 2. We surely have 

Theorem 2. Lei r > 2 and {fi,---,f r } C F and /et Then 
YH=i a rW k n (f s ) converges a.s. to a random variable Y7 s =i a rW^(f s ) 
where each has the same law as W(f s ) found in Theorem.. 

This directly results from Theorem [T] and does not need to be proved. 

3. Application to Extreme value Theory 

3.1. Asymptotic results in the Weibull case. The following sto- 
chastic process 

j=k(n) 

(3.1) T n (f)= /(j)(logX n _ J+1 , n -logX n _ Jin ) 

i=i 
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was introduced by Deme et al.(2012) [3] as a generalization of the Diop 
and Lo continuous generalization of the Hill statistic for f(J) = j T , r > 
1. The Hill statistic corresponding to r = 1 and will be denoted here 
by T n (l). This latter plays a key role in Univariate Extreme Value 
Theorem (UEVT). 

This theory has its foundations in finding the asymptotic law of 
the maximum observation X n<n = max(Xi, ...,X n ). It is said that the 
underlying distribution function F of the observations is attracted to 
some df H if for some sequences (a n > 0) n >i and (6 n ) n >i, we have for 
any continuity point of H, 

lim P ( Xn > n ~ bn < x )= lim F n (a n x + b n ) = H(x). 

It is known that, when it is nondegenerated, H can be parametrized 
as Gy(x) = exp(— (1 + 7x)~ 1 / 7 ),l + ^x > 0, 7 G E named as the 
Generalized Pareto Distribution (GDP). It is said that F is in the 
domain of attraction of G y , hereby denoted as : F G D(G 7 ). The 
reader is referred to de Haan and Feirreira [9] , Resnick [H] , Galambos 
[8] and Beirlant et al. [lj for a modern account of UEVT. 

Although the parameter 7 of the GDP is continuous, the three cases 
(7 < 0), 7 = and 7 > 0, respectively denamed Weibull, Gumbel and 
Frechet cases, may behave ridically differently. But in all the cases, the 
Hill statistic is used to estimate what is called the extremal index in the 
following sense : For 7 > 0, k~ l T n (l) — > 7 as n — > +00 and k/n — > 0; 
for 7 < 0, then the upper enpoint of G(x) = F(e x ) defined by y = 
logsup{a; G R,F(x) < 1}, is finite and A;~ 1 T n (l) /(x -G~ 1 (l-k/n)) -> 
(1 — 7)" 1 n — > +00 and k/n — > and G^ 1 stands for the generalized 
inverse function of G. 

The Diop and Lo generalization of Hill 

j=k(n) 

D n {r) = ^2 ? ( lo S^n-j+i,n - log X n _ j:n ) , r > 0, 

has been studied in ([H]) where its asymptotic normality was proved 
for any 7 but for r > 1/2. Recently, the functional form T n (f), which 
generalizes D n (r) = T n (f T ) = D n (r), has been extendely studied in the 
Frechet and Gumbel cases by Deme et al. ([3]) who proved this : T n (f) 
has a Gaussian limiting process when A(2, f) = Yuj=i fU) 2 /j 2 = +°° 
and 

k 

B n (f) = max{/(j) 2 /j 2 , 1 < j < k}/(J2fU) 2 /f) 1/2 
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— > as n — > oo. It has a non Gaussian limiting process when 
A(2,/)<+oo. 

When particularized for f T , we get that asymptotic normality holds for 
r > 1/2 and not for < r < 1/2. Their results are based on sums of 
independent random variables, and then on Kolmogorov type theorems 
(see [ID])- When put together, for the class of functions f T , we remark 
that the beahavior of T n (f T ) is known for any 7 in the whole extremal 
domain except for the Weibull domain and for < r < 1/2, that if for 
small parameters r's. 

This problem remained unsolved, may be by the fact that its depends 
of sums of dependent data and that we didnot have the appropriate 
frame. Its is clear know that the demimartingale may help to solve it. 
In this paper, let consider the very simple case of 

x - G _1 (l -u)= u lh , < u < 1. 

We use here the index —7 < intead of 7 < 0. We remark that 
G(x) = F(e x ) e D(G_ lh ) if and only if F e D(G_ 1/7 ). We use the 
classical representation of the Yj = logXj associated with the dis- 
tribution function G(x) = F(e x ) through a sequence of independant 
standard uniform random variables Ui, U2; ... , that is 

{Yj,3>\}= d {G _1 (l — Uj),j > 1} 

and then 

{{F n _ j+1 „ n ,l < j < n},n > 1} = d {{G-\l - U j>n ), 1 < j < n},n > 1} . 
This gives 

j=k(n) 

T n (f)/(yo - K-fc+i,n) = 5^ ( lo g X «-i+i,n ~ logX n _ iin ) 

3=1 

= J y^ } f(j) — ~ X ° gXn ~ j ^ ~ ( y ° ~ l °S X n-j+l,n)) 
~[ (Vo ~ Yn-k+l,n) 

j=k(n) 

= d J2 fti) u k,n((u j+ i,n/u k , n y - (u j>n /u k>n y). 

We have for 1 < j < k — 1, 

fc-i fe-i 



(U j)n /Uk, n y = JJiUh, n /Uh+i, n y = exp(-7 r k>g(Uh+l,n/Uh,n 



h=j h=j 
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fe-1 

= exp(- 7 5>fA)- 

h=j 

By the Malmquist representation (see ([15J), p. 336), the random vari- 
ables E^\l <h< n, are independent and standard ones. We arrive 
at 

T n (f)/(y - Y n 

— fc+l,n ) 

k(n) ( fe-l fe-i 

= E /tf) ex p(-^ E E ^i h ) - ex p(-^ E E h ] i h ) 

j=l l ft=j+2 fc=j 



and by (2.1) 



D k , n (f) - T n (f)/(y - K„_fe +1>n ) = (/) 



where Dfe jTl = /(& — 1) — Y^=i fU) s *jk- At this step, we apply Theorem 
[T] to get the final result 

Proposition 1. Let Xi, X2, ... be positive random variable with a finite 
endpoint xo such that log F~ 1 (l — -u) =logxo+M 7 ,0 < u < 1 and 7 > 0. 
Lei / be an increasing function over the positive integers satisfying 
(Kl) and put 

fe-i 

JVC/) = /(*-i)-£ 
T n *(/) = Dfe jn - (r n {f)/(iogxo - io g x n _ fc+1 , n ; 

converges in distribution to a finite random variable W(f) with finite 
expectation. As well, for < r < 1/2,T*(/ T ) converges in distribution 
to a finite random variable W(f T ) with finite expectation. 

We indeed remark that for this simple case in the Weibull case, the 
law of the functional Hill process is determined for < r < 1/2. For 
the general case, we have the following Karamata representation when 
F is in the Weibull case of parameter 7 > : x (F) < 00 and 

(3.2) logxo-F _1 (l-w) = (1 +p(w))w 7 exp( f b{t)r x dt), 

J u 

where (p(u),b(u)) — > (0,0) as u —> 0. In a coming paper, we will 
determine general conditions on /, b and p under which T*(f) bahaves 
as Wj* n as in the present case. 
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3.2. Critical points of the distribution function of W(f). We 
use computer-based methods for approximating the law of W(f). Sim- 
ulation studies show that the empirical distribution functions based of 
BO = 1000 replications are very stable from k = 2000. 

We proceed as follows. Fix r, < r < 1/2, 7 > and k > 2000. At 
each step B from 1 to B0 = 1000, we generate E\{B), E^{B) and 
compute W£ denoted by W%(B). We finally consider the empirical df, 
denoted by Gk, based on W£ (1), W%(B0). Since Gk is stable in the 
sense that it does not significally change from k = 2000, we do approx- 
imate the df Goo of W(f T ) by Gk for k large enough. 

As an example, we illustrate in Figure [T] the df Gk for k = 250, 500, 
750, 1000, 2000, 500 for 7 = 1 and r = 1/4. Here for instance, we fer 
that the support of G^ is [—0.5,0.5]. On the whole, the figures cleary 
establish stability and support our proposal. For users interested to 
use our method, we provide an executable file located at : 



ht t p : / / www . uf r sat . org /lerst ad/resources/ lmhf w 1 . exe 



for the computation of P(W(f T ) < x) = G^x) and P(\W(f T )\ < 
M) = Goodxl) - Goo(- \x\) for x G R. 

3.3. Statistical tests. Let us illustrate here how Goo may be used 
to test the hypothesis that F G D(GLi/ 7 ). We use here the following 
approximation : 



KU) « T nU)l '(log x nin - hgX n _ 



k+l,n) 



We consider gere the statistical test (H) : F G D(G-i/ 7 ) and compute 
the p- values for the models as precised in Table [T} The first three d/'s 
are in the Weibull domain with 7 = 1. The first (Weibull 1) is the one 
we used in the paper. In the two others (Weibull 2 and Weibull 3), we 
introduce a shift of order (1 — u) 5 . For these cases, we conclude that 
the hypothesis is confirmed. But further computations show that for 
n = 300 and k = 200, the hypothesis is rejected for a shift of order 
(1 — u) q with < q < 3. This is conceivable since, as we pointed out 



above, the convergence depends on the functions b and p in (3.2) that 
are here b(u) = and p(u) = u q . This dependence of the results on the 
auxilliary functions will be studied in a coming paper. 



As for the two last cases, the p — values is zero and the related models 
are rejected as expected since the Exponential rv is in the Gumbel 
domain and the Pareto one in the Frechet one. 



10 



ADJA MBARKA FALL, GANE SAMB LO, AND CHEIKHNA HAMALLAH NDIAYE 



k=250 k=SO0 




-0.10 -0.05 0.00 0.05 0.10 0.15 -0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 



k=750 k=1000 




-0.10 -0.05 0.00 0.05 0.10 -0.10 -0.05 0.00 0.05 0.10 



k=2000 k=5000 




-0.05 0.00 0.05 0.10 -0.05 0.00 0.05 



Figure 1. Illustration the distribution functions of Wfe >n (l/4) for 
different values of k 



Models 


Quantiles functions 




T*{fr) 


P-values 


Weibull 1 


F- L (u) = exp(l +u<) 




-0.01897 


50.6% 


Weibull 1 


F- L (u) = exp(l +u<{ 


1 + (1 -«)")) 


-0.02524 


39.6% 


Weibull 2 


F~\u) = exp(l + u f { 


1 + (1-k) 5 )) 


-0.012556 


66.6% 


Standard Expnential 


F _1 (u) = — log(l — u 


) 


0.5313 


0% 


Pareto 


F- L {u) = 1 + U- 1 




0.9969 


0% 



Table 1. Statistical tests for four models using the law of W{l/A) 
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4. Appendix 

This section is devoted to the computations of the moments of 

k-i 

S^exp^Y.Eh/h) 

h=j 

where E' h s are independent standard exponential random variables, 
and their approximations for large values of j. We begin to give a 
particlular and useful for the the expansion of the logarithm function. 
Fact 1. Let e > be fixed for once. There exists < u such that 

< u < uq, log(l + u) = u + 6(e, u)u 2 , 

where 9(e,u) e [— e — 1/2, e — 1/2] = A(e) = [ai(e), a 2 (e)]. For any 
integer m > 1, let Jo (to) such that Jo(to) > ^/(mu ) so that 

j > Jo(m) =>• log(l + 7 /j) = u + 9jU 2 with 9 5 E A(e). 

In the remainder, we concentrate on the moment computations. 

4.1. Moment estimation. 

4.1.1. Exact values. We have for any integer m > 1, 

k-1 k-1 k-1 

E ((S*) m ) = E(exp(- mi Y, E h/h) = Y[E(- mi E h /h) = \[{l+m 1 /h)- 1 

h=j h=j h=j 

E ((S*D = exp ^- Ml + WO j • 
Now for j > Jo (to), 

(k-1 k-1 
-m 7 ^(l//i)-mV^^ 2 
h=j h=j 

4.1.2. Approximated values for moments. We have by 

k-1 



. . , , .1 1 1 . laife)! to 7 
< | fll ( e )| m 2 7 ( 7 - — - j—- 2 ) < ^ 1 



h=j 

For 

|ai(e)| m 2 7 
Ji(£,m) 
we have 

fe-i 



j k-1 (k-1) 2 ' ' 



j > Ji(e,m) V J (to) exp [ -m^^Oh/h 2 j < e £ 

h=j 
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Next by (4.4) 



fc-1 / . v m-y ( fc— 1 j *| 

exp(-m7 = ( t— r J exp(-m7 < ^ t ~ log((fc - l)/j) V 

h=j \ ' {h=j J 



with 

' fc-i 



exp(-l/j) < exp(-m 7 |j] i - log((A; - < exp(-l/(A;-l)). 

We finally have for j > Ji(e, m) V Jo(m), 

5(l,m,j)5(2,m,j), 

with 

< 5(1, j) = 1 + 0( |ai(£)|m27 ) and 5(2, j) = 1 + O^ 1 ))- 

J 

4.1.3. Approximated values for variances. We have for j > Jq{2) 

(fc-i fc-i 
-27^(1//,) -4 2 T 2 E^( 1 )A 2 

and for j > Jo(l) 

(fc— l fc— i \ ^ 

-7E(1A)-7 2 E^(2)A 2 
h=j h=j 

(fc-1 fc-1 
-2 7 ^(l// i )- 2 7 2 E^(2)// i 2 
/i=j h=j 

Thus 

fc-i 

Far(£*) = exp(27^1//i) 

h=i 

{fc-i fc-i 
exp(-4 2 7 2 £ ^(l)/^ 2 ) - exp(-2 7 2 £(26> fc (l) - ^(2))//^ 2 

Since where a; = 4 2 7 2 6 h {l)/h 2 and y = 2 7 2 Yll) h {2)/h 2 are 
both nonnegative, we have \e x — e y \ < \x — y\ . Thus 

fc-i fc-i 
< exp(-4 2 7 2 £ e h(l)/h 2 ) - exp(-2 7 2 £ ^(2))//i 2 ) 
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fc-1 Q 2 

< 2 7 2 £|20 h (l)- 9 h (2)\ /h 2 < 



2fo l£ 

\Z.V h \L) - O h \Z.)\ I U \ ; 

3 



by (4.5). Hence 

(4.2) Var(S*) = (-^) " V(l,j)V(2,j) 

with 

\V(l,j)\ = 1 + O0- 1 ) and < 7(2,j) < 2 fM^)i. 

J 

4.1.4. Covariance approximate values. Let £ > 1 and consider (Tjj+e = 
cov(S* +i , S*). We have 

fc-i 

=exp(^-log(l + 7 /^)) 

j+^-i fc-i 
= exp(^ -log(l + 7//i))exp( - log(l + 7//1)) 

= E(S* +e )eMj2 -log(l + 7A))- 

Also 

fe-i fe-i 

5;5; +£ = exp(-7^^/^)exp(- 7 £ VO 

h=j h=j+i 

j+i-i fc-i fe-i 

= exp(- 7 ^ ^//i-7 Eh/h)exp{-y Yl E ^ h ) 
h=j h=j+e h=j+e 

fc-i j+^-i j+^-i 

= exp(-2 7 £ E h /h)exp{-i ^ = (S* + ,) 2 exp(- 7 ^ E h /h). 

h=j+£ h=j h=j 

Hence 

j+e-i 

E(S*S* +i ) = E (S* +e ) 2 exp( -Ml + 7A))- 
For j > J (l) V J (2), 

«w(sj, = var(s; +£ ) exp ( ]T - i og (i + 1 /h)) 

h=j 
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j+e-1 j+t-l 

cov(S*,S* +i ) = Var(S* +e )eM- 1 £ 1/h - 7 2 £ 9 h /h 2 ) 

h=j h=j 

(4.3) =Var(S* +e ) (j^f^l) " ^ + °(r 1 ))(l+) 

4.2. Integral computations. Let 6 > 1, we get by comparing the 
area under the curve x i — )■ x _b from j to k — 1 and those of the 
rectangles based on the intervals [h, h + 1] , j = 1, .., k — 2, we get 

fc-2 fc-1 

^ h- b < / aT 6 da; < ^ /r 6 , 

ft=j "'J h=2 

that is 

-b 



h=i J 3 



x~ b dx - (k - iy. 



For b = 1, on has 

fc-i 



( 4 - 4 ) r^Y < M(* - !)/i) -£r)< 7- 

For 6 = 2, on has 

(4-5) --^-^E^ 1 
7 fc - 1 7 2 ^ 7 



j fc - 1 j 2 - -j k-1 (k - l) 2 ' 



that is 



1 < In _ j ^ _ STh- 2 < — 
(k-iy- 3 { k-i' ~f 
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